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ABSTRACT 

We find the angular diameter of R Doradus to be 57 ± 5 mas, exceeding that of Betelgeuse and 
implying that R Dor is larger in apparent size than every star except the Sun. R Dor is shown to 
be closely related to the Mira variables. We estimate an effective temperature of 2740 ± 190 K, 
a distance of 61 ± 7 pc, a luminosity of 6500 ± 1400 Lq and a radius of 370 ± 50 R@. The 
characteristics of R Dor are consistent with it being near the edge of a Mira instability strip. We 
detect non-zero closure phases from R Dor, indicating an asymmetric brightness distribution. 
We also observed W Hya, a small-amplitude Mira, for which we find an angular diameter of 
44 ± 4 mas. 

Key words: techniques: interferometry - stars: imaging - stars: individual: R Dor - stars: 
individual: a Ori- stars: individual: W Hya - stars: AGB and post-AGB 


1 INTRODUCTION 

Several nearby red giant and supergiant stars have been resolved 
with 4-m telescopes, producing a number of significant results. 
Angular diameter measurements have provided evidence for first- 
overtone pulsation in Mira variables (Tuthill et al. 1994b; Haniff 
et al. 1995) and produced a new log Tiff versus (J — K) relation 
(Feast 1996). Also important are the discoveries of hotspots on 
a Ori (Betelgeuse) by (Buscher et al. 1990) and of asymmetries in 
the atmosphere of o Cet and similar stars (Wilson et al. 1992; Haniff 
et al. 1992; Tuthill et al. 1994a). 

It has long been assumed that a Ori, a red supergiant, has 
the largest angular diameter of any star in the night sky. Here we 
report measurements of the red giant R Doradus (HR 1492), a 
semiregular variable of spectral type M8, which show it to exceed 
a Ori in angular size. Although the V magnitude of R Dor is only 
5.4, at infrared wavelengths it rivals a Ori as the brightest star 
in the sky. This has lead Wing (1971) to predict a large angular 
size but R Dor has not previously been observed at high angular 
resolution, presumably due to its southerly declination. We also 
report observations of W Hydrae (HD 120285), which is a small- 
amplitude Mira variable that also has a large angular diameter. 


2 OBSERVATIONS 

We have carried out observations of R Dor in the near infrared 
(1.25 fim) and the red (855 nm) using aperture-masking. This tech¬ 
nique involves modifying the telescope pupil with a mask, the ra¬ 
tionale for which has been discussed elsewhere (Haniff & Buscher 
1992; Buscher & Haniff 1993; Bedding et al. 1993; Haniff 1994). 
In brief, an aperture mask with a small number of non-redundantly 
spaced holes provides more accurately calibrated measurements 
at the full diffraction limit. This comes at the expense of lower 
sensitivity and poorer spatial-frequency coverage, and a good com¬ 
promise in the photon-rich infrared regime is to use an annular mask 
(Haniff & Buscher 1992). This has full spatial frequency coverage 
while being minimally redundant: roughly speaking, each baseline 
is measured twice. Furthermore, a thin annulus whose width is com¬ 
parable to the atmospheric coherence length (ro) largely retains the 
advantages of accurate calibration and enhanced resolution. To date, 
however, results from annular masks have only been published for 
binary stars (Haniff et al. 1989). 

In the visible regime, where photon noise generally dominates 
over atmospheric noise, an annular mask is not a good choice. 
Instead, improved sensitivity and spatial-frequency coverage is best 
achieved by using a long slit (Aime & Roddier 1977; Buscher & 
Haniff 1993). 
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Figure 1. Visibilities of R Dor measured at 1.25 (im with the NTT in 1993 
August. The upper panel is for the 7-hole mask, with crosses showing 
visibilities on 21 baselines. The lower panel is for the annular mask, with 
data points being azimuthal averages. In both panels the solid curves show 
the theoretical visibilities expected from uniform disks with diameters of 
52, 57 and 62 mas. 

2.1 NTT Observations 

Observations with the 3.5-m New Technology Telescope (NTT) 
were made with the SHARP infrared camera (Hofmann et al. 1992). 
Details of the experimental setup are given in Bedding et al. (1993, 
1995). We observed R Dor on 1993 August 7 using an aperture mask 
with seven holes arranged in a non-redundant two-dimensional con¬ 
figuration. Each hole had an effective diameter (as projected onto 
the primary mirror) of 25 cm and lay on a circle of diameter 3.05 m. 
The wavelength was 1.25 /im, selected using a standard J filter. 
For both R Dor and a calibrator star (7 Ret), we obtained 500 short 
( 0.1 s) exposures. 

Fringe visibilities were extracted from the averaged power 
spectrum and were calibrated for atmospheric and instrumental ef¬ 
fects by dividing by visibilities from the calibrator star. The upper 
panel of Figure H shows these calibrated visibilities as a function of 
baseline length. It is clear that R Doradus is resolved, despite the 
fact that the visibility at this wavelength does not fall below 50%, 
even on the longest baselines. There is a scatter in the visibility 
measurements at fixed baselines whose cause is unclear, although it 
may be due to a mismatch in seeing conditions for observations of 
R Dor and its calibrator star. Fitting a uniform-disk model to these 
data gives a diameter of (57 ± 5) mas (milliarcsec). 

Note that we expect the calibrator star 7 Ret (HR 1264; V = 
4.5; M4 III) to be slightly resolved. This star has no measured 
angular diameter, but Ochsenbein & Halbwachs (1982) estimate a 
diameter of 11 mas based on its spectral type. Using the V — K 
calibration of Di Benedetto (1993) gives 7.5 mas. Correcting the 



Figure 2. Visibilities of R Dor at 1.25 (im with an annular mask on the NTT 
in 1995 July, showing the average of 11 observations. The solid curve shows 
the visibility expected from a uniform disk with diameter 58.7 mas. 



Figure 3. Visibilities of W Hya measured at 1.45 /im with an annular mask 
on the NTT in 1993 August, based on the average of 9 observations. The 
scatter at short baselines is due to seeing variations between observations 
of the target and the calibrator stars. The solid curve shows the visibility 
expected from a uniform disk with diameter 44.1 mas. 


visibility curve of R Dor for the non-zero size of the calibrator star 
leads to a slight upwards revision in the diameter by 0.5-1 mas, 
depending on the actual diameter of 7 Ret. 

We also observed R Dor on 1993 August 6 using the J fil¬ 
ter, this time with an annular mask that had an effective outer di¬ 
ameter of 3.3 m and a width of 20 cm (lower panel of Figure [l]). 
Further observations using the same setup were obtained almost 
two years later, on 1995 July 19 (Figure 2j). The former data set 
gave a uniform-disk diameter of (57 ± 5) mas and the latter gave 
(59 ± 3) mas, both in good agreement with the results from the 
7-hole mask. The anomalous feature seen at short baselines in Fig¬ 
ure ^ is due to miscalibration caused by slight differences in seeing 
between observations of the object and its calibrator star (Haniff & 
Buscher 1992). The effect is seen to lesser degrees in all of the plots 
presented here, and the affected data were excluded from the fits to 
the visibility function. 

We observed W Hya with the annular mask on 1993 August 8. 
These observations used a variable filter (CVF) with a resolving 
power (A/AA) of about 50, set to a wavelength of 1.45 /im. Two 
calibrator stars were used: HR 5192 and HR 5287. The result, 
shown in Figure [ 3 ), is a uniform-disk diameter of 44 ± 4 mas. 
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Figure 4. Visibilities of R Dor (three observations, upper panel) and a Ori 
(one observation, lower panel) measured at 855 nm with a slit mask on the 
AAT in 1995 January. For comparison, the dotted curves show the theoretical 
visiblities expected from uniform disks with diameters ranging from 50 mas 
(top) to 80 mas (bottom) in steps of 10 mas. 

2.2 AAT Observations 

On 1995 January 13 we observed R Dor and a Ori at 855 nm with 
a slit mask using the 3.9-m Anglo-Australian Telescope (AAT). We 
used the MAPPIT facility (Masked APerture-Plane Interference 
Telescope; Bedding et al. 1994), which consists of optical elements 
mounted on fixed rails at the coude focus. The system differed 
in two important respects from that described by Bedding et al. 
Firstly, the wavelength-dispersed system was not used. Instead, the 
prism and cylindrical lens were replaced by an interference filter 
(A = 855 nm, AA = 40 nm). Secondly, the detector was a CCD 
with on-chip binning, similar to the system used by Buscher et al. 
(1990). 

For each observation we obtained 10000 short (13 ms) expo¬ 
sures of the target followed by an identical number of the calibrator 
star. The calibrator for R Dor was again 7 Ret, while for a Ori we 
used 7 Ori. The aperture mask was a narrow slit whose effective 
width was 8 cm. The slit was aligned diametrically so that the max¬ 
imum baseline corresponded to the full 3.89 m aperture of the AAT. 
The central 1.51 m was obscured by the secondary mirror, which 
results in a gap in spatial-frequency coverage: baselines having 
lengths between 1.19m and 1.51 m are not sampled. 

For R Dor we obtained two observations with the slit oriented 
at position angle 65° on the sky and a further observation at 145°. 
Figure ^ shows calibrated visiblities for these three observations, 
as well as for a single observation of a Ori obtained at position 
angle 115°. The visibilities of R Dor have been corrected for the 


non-zero angular size of the calibrator star, assumed to be 10 mas. 
This correction is only a few percent, even at the longest baselines. 
Note that fringe visibilities from a Ori are likely to be affected by 
the presence of surface features (Buscher et al. 1990; Wilson et al. 
1992). However, our aim here is to establish that R Dor exceeds 
a Ori in angular size, which is demonstrated by our observations, at 
least at this wavelength. We should note that our measured angular 
diameters will be affected by limb darkening and the presence of 
hotspots and non-circularity, which may be different for the two 
stars. By calculating the bispectrum of each observation, we find 
that R Dor shows non-zero closure phases (Figure [ 5 ]). These imply 
an asymmetrical brightness distribution similar to that previously 
found for the supergiant a Ori and Mira variables such as o Cet 
(Wilson et al. 1992; Haniff et al. 1992; Tuthill et al. 1994a). R Dor 
is the first non-Mira giant from which non-zero closure phases 
have been detected. Note that Di Benedetto & Bonneau (1990) 
found anomalously high visibilities on long baselines for the star 
P And (M0 III) which may indicate surface structure, but those 
measurements did not provide phase information. 

The limited position-angle coverage of our data prevents image 
reconstruction and R Dor is clearly a good candidate for more 
detailed observations, particularly in view of its large angular size. 


3 DISCUSSION 
3.1 Angular diameter 

To obtain the true photospheric angular diameter of an M-star from 
a measured diameter requires two corrections, both of which are 
strongly dependent on wavelength. The first correction is for limb 
darkening, which makes a star appear smaller. The second correction 
is for opacity effects from deep molecular absorption bands, in 
which the star appears larger. This latter effect explains why the 
angular size at 855 nm is larger than at J band: the 855 nm filter 
(with 40 nm bandpass) contains a strong TiO absorption feature 
whose depth is about 50% of the continuum (see Fig. 9 of Bessell 
et al. 1989). 

In the near infrared both effects are greatly reduced and the 
measured diameter is within 1-2% of the ‘true’ stellar diameter 
(Bessell et al. 1989, 1996), where the latter is defined by the point 
at which the Rosseland optical depth is unity. We therefore adopt 
our measurement of R Dor at 1.25 /j.m (57 ± 5 mas) as the true 
angular diameter, to be used in the next section. 

Measurements at visible wavelengths can be corrected for 
limb-darkening and opacity effects, although this correction is 
somewhat ad hoc due to the lack of reliable model atmospheres. 
Haniff et al. (1995) made these corrections to their observations 
of W Hya, which were obtained with a slit mask on the WHT at 
wavelengths of 700 and 710nm. They arrived at a photospheric 
diameter of 46 ± 6 mas (after corrections of 30-40%). This result 
agrees with our 1.45 /jm measurement of44±4mas, giving support 
to the accuracy of their correction process. 

For comparison, the uniform-disk diameter of a Ori at 2.2/rm 
has been measured to be 44 mas (Dyck et al. 1992). Thus, while 
R Dor is the largest star in the sky, W Hya rivals a Ori for second 
place. 


© 1996 RAS, MNRAS 000, 000-000 







4 


T. R. Bedding el al. 

3.2 Bolometric magnitude and effective temperature 

The effective temperature of a star can be determined from its 
angular diameter <f> and its apparent flux (Ridgway et al. 1980; 
Bessell et al. 1989): 

log(T ef f/K) = 4.22 - O.lmboi - 0.5 log(0/mas), (1) 

where mboi is the bolometric magnitude. 

Several infrared magnitude determinations for R Dor are avail¬ 
able in the literature (Neugebauer et al. 1971; Kerschbaum & Hron 
1994), which we have converted to the Carter system (Carter 1990; 
McGregor 1994). These are supplemented by several new obser¬ 
vations made using the 0.75-m telescope and the MKII Infrared 
Photometer at Sutherland Observatory, South Africa. To allow ob¬ 
servation of such a bright star, an aperture mask was placed over 
the telescope for some of the observations. The accuracy of these 
measurements is about 0.03 magnitudes at JHK and 0.05 at L on 
the Carter system 

We derive the following average magnitudes: J = —2.51, 
H = —3.50, K = —3.91 and L = —4.30. The scatter between 
measurements at different epochs is about 0.1 to 0.2 magnitudes, 
presumably reflecting the variability of the star. In addition, we have 
the 12 and 25 micron fluxes from the IRAS Point Source Catalogue, 
which are 5160 and 1590 Jy. From these six wavebands we estimate 
the average bolometric magnitude of R Dor to be —0.89 ±0.1 
using a spline-fitting program kindly provided by Dr. Whitelock. 
This ignores the stellar absorption bands in the near-infrared and 
may therefore overestimate the brightness by up to 0.1 magnitude. 

We also observed R Dor in two bands at Sutherland on 1993 
August 7, contemporaneously with our infrared diameter measure¬ 
ment. We obtained J = —2.59 and L = —4.37, implying that 
R Dor was about 0.07 magnitudes brighter than the average values 
given above. On this basis, we take the bolometric magnitude at the 
epoch of our diameter measurement to be —0.96, which yields an 
effective temperature of2740± 190 K. This is somewhat higher than 
previous indirect estimates: Neugebauer et al. (1971) gave 2400 K, 
de Jager et al. (1988) gave 2365 K (after correcting for an obvious 
misprint) and Judge & Stencel (1991) gave 2230 K. The differences 
may be due to the inadequacies of the indirect methods used to 
estimate temperatures front colour of such red stars, and also to the 
intrinsic variability of R Dor. 

3.3 Classification, distance and radius 

R Doradus is an M8 giant and is the brightest, and presumably 
closest, star with such a late spectral type (Wing 1971). The late 
spectral type makes it likely that R Dor lies on the asymptotic 
giant branch (AGB) rather than the red giant branch. The star is 
catalogued by Kholopov et al. (1988) as a semiregular variable of 
type SRb with a period of 338 days. Despite this classification as 
semiregular, R Dor is in many ways closer to the Miras than to other 
SRb stars. Its period is near the peak of the Mira period distribution 
function (250-350 days), while SRb stars almost always have much 
shorter periods (e.g., Kerschbaum & Hron 1992; Jura & Kleinmann 
1992b). Its late spectral type and large J — K are also more typical 
of Miras than other M-type stars (e.g., Feast 1996). 

Further evidence for the Mira-like nature of R Dor comes 
from evidence for mass loss, which is commonly seen in Miras 
with periods of more than 300 days (Jura & Kleinmann 1992a). 
In R Dor, evidence for mass loss comes front the 12-micron IRAS 
flux, which exceeds the ground-based A'-band flux by a factor of 
two. indicating that there may be extended emission. There may 


also be an extended component at 60 and 100 micron, as claimed 
by Young et al. (1993). In addition, the IRAS LRS spectrum (Volk 
& Cohen 1989) appears to show a weak silicate feature, indicating 
recent mass loss. 

Assuming that R Dor is closely related to the Mira variables, 
as seems likely from the preceding discussion, we can apply the 
period-luminosity relations for Miras in the LMC given by Feast 
(1996). For this we use the average magnitudes obtained above. We 
obtain a luminosity for R Dor of 6500 ± 1400 Lq and a distance 
of 61 ± 7 pc, with the result being the same whether we use the 
period-bolometric magnitude or the period-A'-magnitude relation. 
The uncertainty is based on the scatter in the observed period- 
luminosity relation and the uncertainty in the distance to the LMC. 
Our distance for R Dor agrees with estimates of 60 pc by Judge & 
Stencel (1991) and 51 pc by Celis (1995). 

Comparison with the Hipparcos distance will soon be possible. 
Being one of the few red giants accessible to Hipparcos, probably 
closer even than o Cet, R Dor will be valuable in establishing the 
zero-point of the Mira period-luminosity relation. We note, how¬ 
ever, that the accuracy of Hipparcos parallaxes for these stars may 
be compromised by time-varying surface features which, for the 
closest objects, may affect the centroid at a level of a few mas. 

This distance together with our observed angular diameter im¬ 
plies a stellar radius of 370 ± 50 Rq. From the pulsation equation 
Q = P(M/MqY^ 2 (R/Rq)~ 3 ^ 2 and assuming Q = 0.04days 
(appropriate for first overtone pulsation), we derive a mass of 
0.7 ± 0.3 Mg. All the derived parameters are consistent with a 
classification of this star as Mira-like, with the effective tempera¬ 
ture being slightly higher than the average for Miras (Haniff et al. 
1995; Feast 1996). 

Although angular diameter measurements of Miras have 
favoured first overtone pulsation (Tuthill et al. 1994b; Haniff et al. 
1995), there is recent evidence from cluster long-period variables 
in the LMC that the dominant mode is the fundamental (Wood & 
Sebo 1996). If we assume that R Dor pulsates in the fundamental 
mode, which means adopting a value of Q = 0.09 days (Fox & 
Wood 1982), we obtain a mass of 3.5 ± 1 Mq. This is higher than 
expected from its period (Feast 1989; Vassiliadis & Wood 1993), 
giving further support to first-overtone pulsation. 

All previous measurements of the radii of Miras fall in the range 
400-500 Rq, which is taken by Haniff et al. (1995) as evidence that 
Miras are associated with a well-defined instability strip. The fact 
that R Dor shows a more irregular pulsation behaviour but with 
many characteristics of a Mira is consistent with it lying near the 
edge of such a strip. 
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Figure 5. Bispectra from observations of R Dor (lower) and the calibrator 7 Ret (upper) using a slit mask on the AAT. The amplitude and phase of the complex 
bispectrum at each point are represented by intensity and colour, respectively, with blue corresponding to zero closure phase. The axes are spatial frequencies 
along the slit ( u\ and U 2 ), so that each point in the plane defines a triangle of spatial frequencies (in, uo and u\ + U 2 ). Only the nondegenerate portion of the 
bispectrum is shown here, which is contained in a triangle bounded by the lines u\ = 0 , u\ = U 2 and u\ + U 2 = u ma \, where Umax is the maximum spatial 
frequency sampled by the slit (see Buscher & Haniff 1993). The green region in the bispectrum of R Dor at long baselines (high spatial frequencies) indicates 
non-zero closure phases. 


© 1996 RAS, MNRAS 000, 000-000 











